PROBABILITY OF PHASE SEPARATION FOR THE 
BOSE GAS WITH DELTA INTERACTION 



ALEXANDER R. ITS 

Department of Mathematical Sciences, Indiana University-Purdue University at Indianapolis (IUPUI) 

Indianapolis, IN 46202-3216, U.S.A. 
E-mail: itsa@math.iupui.edu 

and 

VLADIMIR E. KOREPIN 

Institute for Theoretical Physics, State University of New York at Stony Brook 
Stony Brook, NY 11794-3840, U.S.A. 
E-mail: korepin@max. physics . sunysb . edu 

and 

ANDREW K. WALDRON 

Institute for Theoretical Physics, State University of New York at Stony Brook 
Stony Brook, NY 11794-3840, U.S.A. 
E-mail: wally@max.physics.sunysb.edu 



ABSTRACT 

We consider the quantum Non-linear Schrodinger equation idt^ = —d£& + 
2c<3>tvl> 2 with positive coupling constant c varying from zero to infinity. We 
study quantum correlation functions of this model using the determinant repre- 
sentation of these correlation functions. We consider the case of a finite density 
ground state and evaluate, using the Riemann Hilbert problem, the asymptotics 
of the probability that no particles are present in the space interval [0, x] in the 
large x limit. We call this the probability of phase separation or alternately the 
emptiness formation probability. 



1. Introduction 



The one-dimensional Bose Gas is described by quantum Bose fields ty(x,t) with 
canonical equal time commutation relations $?(x),W(y) = 5(x — y), The field $ 
annihilates the Fock vacuum \1/|0) = and the Hamiltonian of the model is 



H = J dx [d x m ] d x m + c^ 1 "^} 



The equation of motion is the celebrated nonlinear Schrodinger equation id t ^> = 
—d^ty + 2c\E^\I /2 . The Hamiltonian and the number operator Q = J dx^^ commute 



whereby the model for a given number of particles N is equivalent to the quantum 
mechanical Hamiltonian with delta interactions 



N Q2 

^ = -E^ + 2c E s( Xj - Xh ). (2) 

j=l ux j l<j<k<N 

The Bethe Ansatz for this model was constructed by Lieb and Liniger 0. 

The ground state for the system in a box of length L containing N particles in 
the thermodynamic limit (N — > oo, L — > oo, N/L = D = const.) is a Fermi "sphere" 
of radius q. Introducing the density of particles in momentum space p(A) so that 
D = j; = Jj[ dXp(X), the ground state at zero temperature is then described by the 
linear integral equation 

27 = P(A) ~ 27 F- q ^ (A ' = P(A) ~ 27 (^ P ) (A) ' (3) 

where the kernel of the integral operator K is 

Ar(A '" ) = c» + (A -,.)» • W 

Let us now lay out the structure of this paper. We begin in section 2 by defining 
the probability of phase separation correlation function P(x) which we express in 
terms of the determinant of a particular Fredholm integral operator. In section 3 we 
use the identity ic/(\ — p + ic) = c / °° dse™^^ 1 ^ to write this integral operator as a 
special Fredholm integral operator related to a Riemann-Hilbert problem (RP). The 
introduction of the s-integration means that this RP is infinite-dimensional, or more 
precisely, integral operator-valued. In section 4 we present a new representation of 
Gl(2, C) in terms of such integral operators which reduces the jump condition ([24]) 
of the operator RP to the jump condition (f46|) of a 2 x 2 matrix- valued RP. Finally 
in section 5 we solve the 2x2 RP for large x. Indeed, we find that this 2x2 RP is 
related to that encountered if one studies the case c = oo for which the s-integration 



is unnecessary. The results obtained in section 5 lead us to a conjecture (see [16]) for 
the asymptotic behaviour of the phase separation probability correlation function. 
Presently we are unable to make this conjecture completely rigorous. 



2. Determinant Representation for Probability of Phase Separation 

Although the ground state is homogeneous and translationally invariant, due to 
quantum fluctuations there is a non-vanishing probability that no particles will be 
found in the space interval [0, x] which we denote 

P(x) = Probability that no particles are found in the interval [0, x]. (5) 



In the thermodynamic limit the probability of an empty interval may be expressed 



as the ratio of two Fredholm determinants 



= (0|det(/ + V^ 
det(/ - K/2ir) 

The integral operator K appearing in the denominator is given above in (|j) and is 
x independent. Therefore when calculating the large x asymptotics we need consider 
only the numerator (0|det(/ + V")|0). The definition of the integral operator V is 
more involved. V acts on the interval C = [—q, q] in the same fashion as K, namely 



(Vf)(X) = r dfiViX^)/^). (7) 



However in addition to depending on the parameters x and c (the length of the empty 
interval and coupling respectively) V and in turn det(J + V) are functionals of a dual 

O A 

quantum field V 7 (A). Let us give the kernel of V 



V{\n) = -f 

Z7T 



exp{ix\/2+ ¥ (A)/2}exp{-ix/i/2- V (fj)/2} 
(A — /i)(A — /i + ic) 

exp{zx/i/2+ ¥ (/i)/2}exp{-zxA/2- V (A)/2} 



(/i — A)(/i — A + ic) 



The dual field acts in a Hilbert space with "dual" Fock vacuum |0) so that (0| det(/ + 
V^)|0) denotes the dual vacuum expectation value of the dual field valued determinant 

^ o 

det(J + V) which may be computed from the definition of the dual field V 7 (A) given 
below: 

£ (A) =p (A)+ q (A), (9) 

where V A, /i G C 



c 2 



l^^^^-^-^j, (10) 

[P (X),P (fi)\ = 0=[q (X),q(fi)] ; P (A)|0) = = (0| q (A). (11) 
From the symmetry of log(c 2 / (c 2 + (A — /i) 2 )) in A and fi we have the crucial property 

(A), V (//)] = VA,/iea (12) 



This allows us to treat <P (A) simply as some function on C in our calculation of 
det(I+V). Only at the end of this calculation do we take the dual vacuum expectation 
value. 



Obtained using the Riemann Hilbert problem in the large x limit, our conjecture 

o 

for the determinant as a functional of V 5 (A) is a : 



det (I + V) *=°° (const) exp { -(xq) 2 /8 - ± F + o(x) 1 . (13) 

2tt J-q ^/j, 2 - 



Let us now analyse the dual vacuum expectation value of (fl"3|). Using (P[)-([TTD and 
the identity e A+B = e A e B e -^ A ' B ^ valid V A, B which commute with [A, B] we obtain 

, . , r , , f (xq) 2 x 2 ft ada ri vdv , c 2 + (u — z/) 2 ] 

det J + V) 0) ~ exp J_i«L _ / ———log- 1 



87T 2 7- g ^/i 2 - 9 2 7-9 - 5 2 C 



.2 



2 



expj-^[l + /(c 2 /g 2 )]j, (14) 

where we denote the dimensionless integral 

r / 2/2 ^ 2 f 1 ydy f 1 zdz f c 2 /q 2 + {y + z) 2 \ 

/(C /q ) = ~ 2 L 7T^¥ L bg { c 2 /q 2 + (y- Z ) 2 ) • (15) 

I(c 2 /q 2 ) has the following properties: I(c 2 /q 2 ) > 0, g^p) < 0, I(c 2 /q 2 — > oo) = 
and I(c 2 /q 2 — > 0) = 1. Therefore we find as our conjecture for the phase separation 
probability correlation function 

"=°° (const) exp | --^2 [1 + /(c 2 /g 2 )] +o(x 2 )| . (16) 

The probablility of phase separation is Gaussian in the length x and the rate of decay 
in x decreases monotonically as the coupling c varies from zero to infinity. 



3. Integrable Linear Integral Operators and s-Integration 

Integral operators I + V (acting on some interval C say) where the kernel of V 
has the form 

1 N 

V(\,») = - 5>(*)i5i(A0, (17) 

A - li j=1 

with = ^2j =1 ej(X)Ej(\), belong to a special class which we will call "integrable" 
operators (the functions ej(A) are N linearly independent functions, continuous and 
integrable on the interval C, similarly for Ej). Indeed such integrable integral op- 
erators form a group (where multiplication is defined in the usual way for integral 

a + ^/ A 2 — q 2 denotes the limit approaching the contour C from above. Further, we remind the reader 
that the Hardy small "o" symbol o(x) denotes any function of x decreasing faster than x as x — > oo. 



operators) and are closely related to a certain N x N matrix Riemann-Hilbert prob- 
lem. 

Therefore to employ the RP for finite c we must express the kernel (El) in the 
standard form above. To this end we rewrite the second factor in the denominators 
using an old identity made famous by Feynman x _^ +ic = —i Jo° dse ls ( x ~ fJ " + ' iC \ The 
kernel (§) becomes 

V(X,fi) = r e4Ms)e-^\s) ~ e.(Ms)eMs) dg (lg) 

JO A — fl 

°° e 1 (X\s)E 1 ( f i\s) + e 2 {X\ s)E 2 {ii\s) ^ 



o A — ii 

where e±(A|s) = J— exp{±(ixA/2+ ¥ (A)/2)}v^exp{±isA - cs/2}, (20) 
V 2n 

and (e a (A|s),e 2 (A|s)) = (e+(A|s), e_(A|s)), (21) 

(E 1 (X\s),E 2 (X\s)) = (e_(A|s),-e + (A|s)), (22) 

which is in the form (|T^) with the sum J2f=i generalized to an integral f£° dsY%=i- 
We must, therefore study a 2 x 2 operator valued RP which we now define. 

The operator valued RP is to find the 2x2 matrix x(A) whose entries are integral 
operators acting on the interval [0, oo) and functions of the complex plane C with 
kernels 

X(AM= , XeC;s,te [0,oo), (23) 

V X2i(A|s,t) X22(X\s,t) I 

such that: 

1. x(A|s, t) is analytic V A G C \ C. 

2. x(A|s,t) is, however, not analytic across the conjugation contour C, rather the 
limits x+(A|s,t) and x_(A|s,t), approaching C from above and below respec- 
tively^ satisfy the jump condition 

x ^(X\s,t)= J2xi(Ms,r)G jk (X\r,t)dr, V A G C, (i, k = 1, 2), (24) 

where the conjugation matrix Gik(X\s,t) is given by 

G ik (X\s,t) = 5 ik 6(s-t)-27rie i (X\s)E k (X\t), (i, = 1,2). (25) 



fc We consider C = [—q,q] as an oriented contour so that J c = f_. In general for any oriented 
contour we will define "+" to be the limit from the left side of the direction of travel. 



3. x{M s >t) 1S canonically normalized: 

X(°°M= f J J (26) 
In fact, if x(<\) solves the RP defined above, then at A = oo it has the whole 



1 



Laurent series expansion x(M s > t) = ^ g ^ J $( s ~ t) + E™=i Mn ^'^ where the 

integral operators M n are independent of A. 

In order to extract the determinant det(J + V) from the operator RP we consider 
the formal solution to the above RP 

(\\ A A s( .x , f fi(fJ.\s)E k ((i\t)d(j, 

Xifc(As,t) = d ik d(s-t)+ (27) 

Jc /i — A 

where e;(A|s) = [(/ + t>)/,](A|s) = /,(A| S ) + / V(X, n)fi(ii\s)dii. (28) 
For A large we have 

Xtfc(A|s,t) A ~°° ^5( S -t)-^| c /,( /U | S ) J E; fc ( / i|t)rf/i + 0(A- 2 ) (29) 
= ^-t) + M ^ g,t) +0(A- 2 ). (30) 

But now (see i) calculate the logarithmic derivative w.r.t. x of det(J + V) from (19) 
d "if f°° 

— logdet(/ + V) = - (fi(ii\s)ez(ji\s) + h(fi\s)ei(ji\s))dsdii 



dx 2 Jc Jo 



2 Jo 



(M 1 u (s, s) - M 22 (s, s))ds = -itiMl 1 , (31) 



where "tr" denotes the trace operation in the space of integral operators on [0, oo). 
Therefore det(/ + V) can be obtained from the large A asymptotics of the solution to 
the above RP. 

At this point the prospect of solving an operator valued RP should seem daunting 
since even the solution for the matrix case is not in general known. However we now 
construct a representation of Gl(2, C) which will allow us to consider a (known) 2x2 
RP. 

4. Gl(2, C) Represention by Integral Operators 

Let us consider some 2x2 matrix of integral operators O with kernels 

/ e>n M) o 12 ( s ,t) 

V 021 M) o 22 { s ,t) 



0(s,t) 



s,t e[0,oo), (32) 



and multiplication denned by c 

O m (s,t) = (OiOn) (s,t) 



0$(8,t) 0™(s,t) 
0™(s,t) 0™(s,t) 



s,t e [o, oo) 



(33) 



where 



On(*,t) = J~dr (O I 2 i(a,r)0? 1 (r,f) + Ol 2 (s,r)Og(r,t)) 

roo . 



(34) 



Let us now construct a special class of such operators O which form a representation 
of GZ(2,C) (the generalization to Gl(N,C) is straightforward). Consider a pair of 
functions ee(s) and fi(s) on [0, oo) which we write in Dirac notation as a(s) = |1) and 
(3(s) = 1 2). In this notation we may write left multiplication by Oik as 



O, 



re 

1 = / 

Jo 



dtO ik (s,t)a(t). 



(35) 



Further suppose the functions A(s) = (1| and B(s) = (2| (right multiplication by Oik 
is defined integrating over the first argument of Oik(s,t)) satisfy 



roo roo 

(1|1) = / dsA^a^s) = 1 = (2|2) = / dsB(s)P(s). 
Jo Jo 

Observe now that one may define a representation A of Gl(2, C) via 



(36) 



Me GZ(2,C) i — >■ A(M) = 







V 



o 



/ M n |l)(l| M 12 |l)(2| \ 



J-|2)(2| J 



\ M 21 |2)(l| M 22 |2)(2| J 

(37) 

(Mn, Mi 2 , M 2 i and M 22 are complex numbers and / is the identity operator in 
the space of integral oprators on [0, oo)). Multiplication by the integral operators 
(projectors) |1)(1|, |1)(2|, |2)(1| and |2)(2| is given, for example, by 



ii><2i/( a ) = (jr 00 d S s( a )/( a ))ii>. 



(38) 



c We assume that the integral operators Oik act on the interval [0,oo) but this is of course not 
necessary. 



In particular, we have [I - |1)(1|] |1)(1| = 0. Indeed for any M, N e Gl(2,C) the 
representation A has the following properties: 

A(MN) = A(M)A(N) ; A(I) = I ; A^M^ 1 ) = A~\M) (39) 
det A(M) = det M = M n M 22 - M 12 M 21 (40) 



tr A(M) - 








trM = M u + M 22 (41) 



Notice that ( |4"0"D expresses an infinite dimensional determinant as the determinant of 
the 2x2 matrix M. Let us conclude this section by rewriting the conjugation matrix 
of our RP using the representation A. We begin by defining functions 



|1) = ^ e lsX - cs > 2 = a(s), (1| = ^e~ lsX - cs / 2 = A(s) 
|2) = ^e~ isX - cs ' 2 = P(s), (2 1 = ^e isX ~ cs ' 2 = B(s) 



(42) 



The vectors |1) and |2) are normalized but not orthogonal, so that their inner products 
are = cf™e- cs ds = 1 = (2|2) and (1|2) = (2|1)* = c/ °° e - 2isX - cs ds = 

(for A real the bra and ket notation denote complex conjugation in the usual fashion 
although this is not necessary). Their outer products are of course integral operators 
in terms of which we will write the conjugation matrix. We now have 



ei(A|s) = J± e i™W+nW)\i) = e °(A)|l>, 

E^t) = ^T e (-teA/2-V(A)/2)^| = £°(A)(1|, 

e 2 (A| S ) = v /i e C--A/ 2 -^(A)/ 2 )| 2> _ e o (A) | 2)j 
^ E 2 (X\t) = -i/Je(« A /wW/ 2 )(2| = £ 2 °(A)(2|. 



(43) 



Hence, from (|25| ) we have for the conjugation matrix (where from now on the s and 
t dependence will be taken as understood) 



G(X) 



I -27rie° 1 (\)El(\)\l}(l\ 
-27rze°(A)E 1 °(A)|2)(l| 



-27u e °(A)£°(A)|l>(2| 
27ue°(A)£ 2 °(A)|2>(2| 



A(G°(X)). (44) 



where the 2x2 matrix G°(X) which here plays the role of M is readily computed to 
be 



G°(X) 



( 2 _ e ix\+V{\) \ 

t e -ix\-!p(\) Q j 



(45) 



Except for the additional exponents ± *P (A), G°(X) is the conjugation matrix one 
finds for the 2x2 RP encountered in the case c = oo. 



5. Asymptotic Analysis of the Riemann Hilbert Problem 

Firstly let us sketch the derivation of the asymptotic solution to the 2x2 RP 
with conjugation matrix G°(X) given above in ( jig) (we follow the scheme developed 
in i for the case V 9 (A) = 0). I.e. we are searching for the 2x2 matrix function x°(A) 

normalized to ^ ^ J j at infinity and analytic everywhere in the complex plane C 

save for a discontinuity across the contour C = [—q,q] where x°W satisfies 

X°-W=X+WG°W, VAeC. (46) 
The idea is to make a "gauge transformation" under which the above RP limits for 
large x to a new RP with conjugation matrix G = ( ^ 0^ ) ' ^ e ^ ne ' t nere f° re 

X°(A) = exp {^3} x°(A) exp {[(ix/2)(X - g(X)) - ^(A)]a 3 } , (47) 

(03 = diag(l, — 1)). The function g(X) = a/A 2 — q' 2 satisfies the scalar RP <7_(A) = 
—g + (X) across C with normalization g(X) A ~°° A — q 2 /2X + 0(A~ 3 ) and 



2m Jc ^/fx 2 — q 2 (fi — X) 

O 

satisfies the affine scalar RP + ip-(X)+ V 9 (A) = across the contour C with 

o 

normalization ip(X) — ipoo = ir~- Jc lf ")^ d ^' . Notice that ip(X) is well defined on all 

o 

of C even though it depends on the dual quantum field V 9 (A). One may verify that 
x(X) satisfies a canonically normalized RP across C with conjugation matrix 

G (A)=exp{-[(^/2)(A-^ + (A))-^ + (A)]a3}G (A)exp{[(^/2)(A-^4A))-^(A)]a 3 } 



/ 2e^ i - X)+ix9+i ' X) -1 



x— >oo 



/ -1 



\1 



(49) 



where Aifj(X) = — ip-(X). The canonically normalized RP for X(X) is easily 

solved in the large x limit, we find 



X°(A) 



/ |(a(A) + a- 1 (A)) |(a(A) - a'^X)) 
v -|(a(A) - a- x (A)) |(a(A) + a" 1 (A)) 



(50) 



where a(A) = (j^) 1 ^ satisfies the canonically normalized RP a_(A) = ia + (X) (VA G 
C). Inverting the transformation ( |47| ) we obtain the large x solution to the original 
RP, 



X°(A) 



-i/(AH' 00 -Hja;(A-g(A))/2 * 



/ i(a(A)+a-\A))e^ A ^^ A ^ A »/ 2 f(a(A) -cT^e- 
^-f^A) -a-\A))e^ A ^^ A -^ A »/ 2 ^a(X) + a~\X))e-^ x ^ Mx ^ x ^/ 2 / 



-1 



V 



' 2 C 



—-01 + ixq 2 /4: 



0(A 



-2> 



(51) 



where V>i = ^7 /c 



We may now try to utilise our GZ(2, C) representation to write the solution to the 



finite c operator valued RP as 



since by property 



X (X)=A( X °W), 



(52) 



X_(A) = Atf_(\)) = ^( X +(A)G°(A)) = A( X ° + W)A(G°(X)) = %+{\)G(\). (53) 

The observant reader may notice that this solution does not obey the normalization 
at infinity since the projectors |1)(1|, |1)(2|, |2)(1| and |2)(2| are essentially singular 
at A = oo. Moreover, as bounded operator valued functions they are analytic in the 
strip |ImA| < c/2 only. Therefore strictly speaking equation ([5^) is not correct. The 
correct version is 

x(A) = HX)A(x°(X)), (54) 

where $(A) is an operator valued function analytic in the strip S : |ImA| < c/2. 
The last equation itself may be treated as a Riemann Hilbert Problem. Indeed, 
let T be a simple contour belonging to the strip S encircling the interval C in the 
counterclockwise direction. Introducing 



Y(X) 



x(X) if A is outside V 
$(A) if A is inside T, 



(55) 



we can rewrite our initial operator RP 1-3 as the following new RP d : 
l r . Y(X) is analytic VA G C \ T. 

2 r . Y-(\) = y + (\)A(x°W), A er. 

d To avoid confusion, in agreement with our conventions let us reiterate that + is the limit from 
inside the anticlockwise contour T. 



3 r . f(oc)=J ( I I 



From ( pip and the equation 

det xn(A) = 1 + tiM^/X + 0(\- 2 ), (56) 

we see that to calculate det(J + V") we need det Y U (X) for large A. It follows from 2 r 
that 

det Y + n (A) = detf n (A) ( X ° 2 (A) - x^(A) C (A)) , (57) 

where "det" indicates the determinant of integral operators on the space [0, oo). and 
C (A) = (l|[fi 1 (A)]- 1 Y_ 12 (A)|2) = (l|[x 11 (A)]- 1 x 12 (A)|2). Let us now assume that 

(a) c(A) is analytic outside T, 

(b) c(A) = 0(1/ X) as A -> oo. 

Note that the product [fi 1 (A)]" 1 r 12 (A) (= [x"(A)] -1 xi 2 (A)) has these properties. 

Statements (a) and (b) constitute our main hypothesis concerning the solution 
x(A) of the RP 1-3. By virtue of this hypothesis, equation ( p7[ ) implies that 

det Y n (A) = ( X ° 2 (A) - X2i(A)c(A))-\ (58) 

for all A outside of the contour T. In addition property (b) leads to the equation 

lim A(detF n (A) - 1) = - lim A( X ° 2 (A) - 1). (59) 

A— >oo A^oo 



Taking (j57|) and ([H]) into account we obtain 

trM* 1 x =°° ij} X - ixq 2 /A + o(l), (60) 
which yields our conjecture for the determinant 

det (/ + V) *=°° (const) exp |-(xg) 2 /8 - f - / + o(x)\ , (61) 

I 2n Jc ^/fi 2 - q 2 j 

o o 

as quoted in section 2. Finally let us observe that for c = oo, [p (A),<? (/z)] = so 

o o 

that for any functional F^P (A)] the dual vacuum expectation value (0|F[y? (A) 10) = 
F |o Hence (0| det(J + V)\0) | c=00 = e - (xg)2/8 . An old result first obtained inB and 
then extended to the complete asymptotic expansion in i. 
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